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Abstract 

Statistical evaluations of the economic mobility of a society are more difficult than measurements of the income 
distribution, because they require to follow the evolution of the individuals’ income for at least one or two generations. 
In micro-to-macro theoretical models of economic exchanges based on kinetic equations, the income distribution 
depends only on the asymptotic equilibrium solutions, while mobility estimates also involve the detailed structure 
of the transition probabilities of the model, and are thus an important tool for assessing its validity. Empirical data 
show a remarkably general negative correlation between economic inequality and mobility, whose explanation is still 
unclear. It is therefore particularly interesting to study this correlation in analytical models. In previous work we 
investigated the behavior of the Gini inequality index in kinetic models in dependence on several parameters which 
define the binary interactions and the taxation and redistribution processes: saving propensity, taxation rates gap, 
tax evasion rate, welfare means-testing etc. Here, we check the correlation of mobility with inequality by analyzing 
the mobility dependence from the same parameters. According to several numerical solutions, the correlation is 
confirmed to be negative. 
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1 Introduction 


The rise of economic inequality has become, in recent times, cause for serious concern and the object of frequent 
debates. Inequality is increasing in several countries and also appears to be accompanied by a diminution of social 
mobility. This combination is especially harmful, because it frustrates the attempts of many individuals to move up the 
income ladder and improve their situation. Empirical surveys confirm the existence of a negative correlation between 
mobility and inequality JTjO. It seems therefore very important to highlight the fundamental causes of this correlation. 

Our aim in this paper is to explore these themes through mathematical methods, by providing a suitable definition 
of mobility and by comparing the possible scenarios compatible with different policies and conditions. In this connec¬ 
tion, we introduce some indicators apt to measure various aspects of mobility within a family of kinetic-type models 
for the description of economic exchanges in a closed society. We investigate the relationship between economic 
inequality and mobility by analyzing income distributions emerging from different fiscal and welfare systems. The 
output of a large number of simulations shows clear evidence of a negative correlation between mobility and equality, 
in agreement with the qualitative features of real world situations. The models also allow to study the simultaneous 
dependence of the Gini index and the mobility parameters on taxation rates and welfare means-testing, and to track 
variations in the mobility of different income classes. In particular, it turns out (the simulations suggest) that there 
are, with reference to variables measuring welfare and taxation, lines on which both the Gini index and the mobility 
parameter are constant. 

Our approach fits in with a complex system perspective, in the sense that we find the emergence of collective 
macro-level patterns out of a whole of interactions occurring at the micro level. The need of considering, when 
dealing with socio-economic questions, a multiplicity of individual interactions rather than a representative agent is 
highlighted and advocated e.g. in 01 SI 13 • For this purpose, agent based computational models have been typically 
developed and applied within the economics community, j6j. On the other side, during the last two decades a number 
of physicists have started to tackle economic and financial questions by using tools of statistical physics, see e.g. 
0130 and the recent review itTOll . The reason for that naturally arises from certain analogies existing between 
markets or populations of individuals who exchange money and particle systems or gases in which molecules collide 
exchanging energy. In particular, as illustrated below, our methods and models allow, thanks to their characteristic 
features, the definition and investigation of a mobility concept. 

The remainder of the paper is organized as follows. Section[2]sets up a framework for the description of monetary 
exchanges in a closed society. In Section [3] a family of models is constructed within this framework. In the fourth 


section the concept of social mobility is shortly discussed and quantitative indicators to measure it are introduced. 
Hence, the correlation between the economic inequality and mobility is explored. Finally, a summary is contained in 
the last section. 


2 A discrete-type framework for the description of monetary exchanges in 
a closed society 


Moving from a stylized description of the complex of monetary exchanges which take place in a closed market society, 
we outline in this section a discrete-type framework suitable for its modelling. This framework was proposed and 
discussed in our previous works on Cana, but we give here a slightly different introduction. The problem at issue 
concerns the formation of the income distribution curve and the economic inequality in a population of constant size. 
If the variable r £ [0, +°°) denotes the individual income and x(r) the density of the population having income r, 
one could say that the total population and the total income are given by f 0 °°x(r)dr and Jq °° rx(r) dr respectively. 
However, it is quite natural and also convenient in view of a numerical treatment of the problem to divide the population 
into a finite number of classes. Toward this goal, we consider the 22 + 1 numbers 0 = ro < r\ < r 2 < ... < r„ with r„ so 
large as to represent a sort of supremum of the possible income of any individual in the given society. For simplicity 
we assume the density of the population whose income belongs to the interval [r,-_ i, r ( ) to be constant and denote it by 
Xj for i= 1, .... ii. The fraction of individuals having income in [r, \. r,) is then 


Xi = / xi dr = Xi(n - r,_ i) 

■> r i-\ 


( 1 ) 


and its income is 


Ri = 



rxjdr =x. 



( 2 ) 


We refer to the class of individuals with income in [r,-_i, r,) as to the f-th class and to its individuals as to /-individuals. 
An immediate consequence of 0 and ([2]) is that R, = Aj r,, where r, = (r,-_ 1 + r,j/2 denotes the average income of 
individuals in the i-th class. 


We assume that pairwise monetary exchanges take place between single individuals. A "microscopic" description 
of this can be as follows: when an k-individual pays to a k-individual an amount of money S (with S < (r i+ i — f,j for 
all i = 1 ,...,«), the h-th class gets slightly poorer, while the k-th class gets slightly richer. This induces a retrocession 
of a portion of the population of the / 2 -th to the (h — l)-th class and an advance of a portion of the k-th population to 
the (k + l)-th class. Easy calculations show that, if the total population and the total income have to remain unchanged 
as happens in the conservative case dealt with here, then the variations of population in the / 2 -th and in the k-th class 
are given by 

2 s 2 S 

A Xi, = -and AX^ = -. (3) 

rh-fh-2 n+i-n-i 

How often monetary exchanges can occur can be postulated through the definition of suitable terms p^.k for h,k = 
1 , ..., 22 , with each phk accounting for the encounter frequency rate of individuals of the / 2 -th and in the k-th class and 
expressing as well the probability that in an encounter between an / 2 -individual and a k-individual, the one who pays is 
the / 2 -individual. The possibility that none of the two pays has to be taken into account as well. Hence, for any h and 
k in { 1 , ..., 22 }, the inequalities 0 < p^k < 1 and p) u k + Pk,h < 1 have to hold true. Apart from that, the values of the 
Pk k should be determined based on phenomenological observations. A specific choice for the values of the /;/, /. will 
be done in the third section of this note. 


Beside direct money exchanges, we also assume the existence of a process of taxation and redistribution. If we 
imagine that the money received form the k-individual is subject to taxation with a tax rate equal to Tj. £ [0,1), sup¬ 
posing (at least as a first case study) that the k-individual pays the due tax, we observe that the entire operation is 
equivalent to one in which the / 2 -individual pays a quantity S, the k-individual earns a quantity .S' (1 — T/J, and the 
quantity S Xk goes to the government. We then imagine, here again postulating an irreproachable situation, that the 
government redistributes the tax revenue (which is just S Xk for the single mentioned interaction) to the entire popula¬ 
tion through the provision of welfare services: education, health assistance and so on. At this step, we observe that 
in fact the government intervention can be ignored, since the mechanism can be equivalently described as follows: 
in correspondence to the mentioned direct interaction between an / 2 -individual and a k-individual, individuals of all 





classes receive an amount of money corresponding to a returned tax portion. In other words, the taxation and redistri¬ 
bution effect is equivalent to the effect of an indirect interaction between the paying /z-individual and all others Q Of 
course, the indirect interactions too contribute to the possibility of class change of individuals. 

In order to formalize the overall mechanism described so far we introduce 

— the direct transition probabilities CL £ [0,1] for i,h,k £ {1,...,«}, satisfying 

T J C hk = l for any h,k £ {1,.(4) 

i—1 

where C‘ hk expresses the probability that an individual of the /z-th class will belong to the z-th class after a direct 
interaction with an individual of the k-th class, 

and 

— the indirect transition variations 7j^ : R" —> R for i,h,k £ {1, ...,zz}, where each 7j^ (X) withX = (Xi £ 

R" is a continuous function, satisfying 

n 

for any h,k £ ,n} and X £R”, (5) 

z—1 

accounting for the indirect interactions and expresses the variation in the z-th class due to an interaction between an 
/z-individual and a ^-individual. 

As will be shown in the following section, the construction of both the parameters C l hk and the functions T [hk] is 
accomplished combining the variation terms ([3j) and the probabilities /?/, £. 

The variation in time of the number of individuals belonging to the n income classes is described by the following 
balance differential equations: 

5 = IEW-lt C ^*+I E r hki X Wk for i £ {1, ...,zz}, (6) 

aT h=lk =1 h=lk= 1 h=\k= 1 

which, in view of Q, can be simplified to take the form: 

= [c\ l k + T{ hk] {X))x h X k -X i Y d Xk for i £ {l,...,zz}. (7) 

QT h=lk= 1 k=\ 

They are nonlinear ordinary differential equations. 


3 An example of model family 


In order to construct a specific model or model family, we start by fixing specific values for the ph.k- We take them as 


Ph,k = min{r hl r k }/4r n , 

with the exception of the terms pjj = ?j/2r n for j = 2,...,zz — 1 , p^i = h/2r„ for lr = 2,...,zz, p n _k = fk/2r n for 
k = 1 1, p\,k = 0 for k = 1 , ...pi and p] m = 0 for h = 1, ...,;z. This choice is suggested by the consideration that 

poorer individuals typically earn and pay smaller amounts than richer ones and do so less frequently. It guarantees a 
degree of heterogeneity among people of different classes, also in connection with their interactions with others. 

1 We point out that actually, for a technical reason, in the following we will assume that no redistribution goes to individuals of the «-th class. 
Indeed, if an individual of the n-th class would receive some money, the possibility would arise for him to advance to a higher class, which is 
impossible. 



Recalling now the formulae pL we assume that the only possibly nonzero elements among the C' hk are: 


C‘+1 ,* 

25(1 -T t ) 

— Pi+1 ,k 5 

n+i - 

C U 

, 25(1 -Ti) 

= 1 Pk,i 

r i+\ — r i-l 

cU, k 

25(1 -Zi-i) 

= Pk,i-1 

r, - n- 2 

We stress that the expression for C‘ t k 

in |8| holds true for i < 


Pi.k - 


r, - r t - 2 


( 8 ) 


expression for Cj k is effectively present only provided i < n — 1 and k > 2, while its third addendum is present only 
provided i > 2 and k<n — 1; the expression for C\_ x k holds true for i > 2 and k > 2. 


As for the T' hk , we first introduce certain coefficients w,- for i = 1..... u, which denote the weights of a differently 
distributed welfare. In the following we will take the rj (and consequently the rf) linear in j: rj = cj for some constant 
c > 0. Then, a conceivable expression for the weights w,, according to which, the wj decrease linearly as functions of 
j, is given by 


2 /. n+l\,_ . 

Wj = r„+i-j +—j r ( J — J 


( 9 ) 


with y € (0,1 /2). The smaller ye (0,1 /2), the larger is the difference w i — w„. Indeed, w i — w n — (r n — h) (1 — 2y). 
The value y = 1/2 also is admitted, but if y = 1/2, then wj has the same value for each j = 1 


We can now express the T [hk\ as 


T{ hk] (X)=U' hk] (X) + V’ hk] (X), 


[**] 


[«]' 


wher^e] 


U[ hk] (X) 


Ph,k^Sx k / Wj-iXj-i _ WjXj \ 

L”= 1 W jXj V n - Vi- 2 r i+ 1 - r,-_ 1 J 


( 10 ) 

( 11 ) 


represents the variation corresponding to the advancement from a class to the subsequent one, due to the benefit of 
taxation anc0 

Sh,i +1 < 5 h,i 


WjXj 

V{ hk] (X) = PKk 2ST k J - 1 


I"=t ^Xj 


n,-n -1 n t -n ~2 


( 12 ) 


with 5/, ^ denoting the Kronecker delta, represents the variation corresponding to the retrocession from a class to the 
preceding one, due to the payment of some tax. 


We notice here that the effective amount of money representing taxes, which is paid in correspondence to a payment 
of 5(1 — T k ) and is then redistributed is S T k (J/’I J WjXj)/(Jj" =l WjXj) instead of S T k . This is in agreement with the 
fact that, due to the bound on the admissible income, individuals of the n-th class do not benefit from the redistribution 
(see Note 1). 


Finally, to fix ideas, we take 5=1, 


rj = 2Sj, 


(13) 


and 


j — 1 , , 

tj — tmin H r ( tmux ^min )? 

n— 1 


(14) 


for j = 1, where T„,/„ and T max respectively denote the minimum and maximum tax rate. Still, the value of y, T m; „ 
and z max have to be fixed. Hence, the equations (0 describe a family of models rather than a single model. 


Finding analytical solutions for these equations is in practice hopeless. But, as will be shown in the next section, 
several relevant facts can be understood through simulations. Before passing to numerical results, we would like to 

2 In (TTJ, h > 1 and the terms involving the index i— 1 [respectively, i+ 1] are effectively present only provided i— 1 > 1 [respectively, i +1 < n]. 
In other words: the 1° term into parentheses on the r.h.s. of CD is present for 2<i<n\ the 2° term into parentheses on the r.h.s. of CD is present 
for 1 < i <n — 1. 

3 In (JT2J, h > 1 and the terms involving the index i— 1 [respectively, i +1 ] are effectively present only provided i — 1 > 1 [respectively, i +1 < ri\. 
In other words: the \° term into parentheses on the r.h.s. of 0 is present for 1 < i < n— 1; the 2° term into parentheses on the r.h.s. of 0 is 
present for 2<i<n. 
















mention the following two properties of Q. For both of them the proof of the corresponding result established in IfTTI 
and referring in fact to a slightly different model can be easily adapted. 

Well posedness of the Cauchy problem. For fixed values of the parameters y, T min and T max , in correspondence 
to any initial condition Xq = (Xqi.... ,Xo„), for which Xo,- > 0 for all i = 1, ...,n and L/Li x 0i = 1, a unique solution 
X(t) = (X\ (;).... ,X n (t)) of |7]) exists, which is defined for all t £ [0, +°°), satisfies X (0) = Xo and also 

n 

Xi(t) > 0 for i = 1,and ^Xj(t) = 1 for all f > 0. (15) 

1=1 


Conservation of the total income u . For fixed values of the parameters y, T„,„, and X max , the scalar function p (X) = 
Y!l=\ hXi, expressing the total (and mean) income, is a first integral for the system (J7J). 


From now on, we report on the outputs of a large number of simulations. To run them, we fixed n = 15. The first 
fact one observes is the following one. 

Uniqueness, for any fixed value of p, of the asymptotic stationary distribution. For fixed values of the parameters 
7 , T min and Xmax, for any fixed value n £ [n, r„], an equilibrium of (|7]) exists, to which all solutions of (J7]), whose initial 
conditions Xo = (X 01 ,... ,Xo„) satisfy Xo, > 0 for all i = £" =1 Xb,- = 1, and E " =1 H x 0i = £i tend asymptotically 

as t —► + 00 . In other words, a one-parameter family of asymptotic stationary distributions exists, the parameter being 
the total income value. 


Other properties can be established from comparisons between different models of the family under consideration. 
In previous work ( fll4l 15]) we investigated the effects on the shape of the asymptotic income distribution curve and 
the economic inequality produced by the adoption of fiscal systems with different spread between the maximum and 
the minimum tax rates and also produced by the introduction of differently weighted welfare measures. Briefly, the 
simulations suggested the following properties. 

Dependence of the asymptotic stationary distribution on the tax rate difference x max — T min- The profile of the 
asymptotic stationary distribution depends on the difference between the maximum and the minimum tax rate: if 
T max ~ ?min is enlarged, while all other data are kept unchanged, then an increase of the fraction of individuals belonging 
to the middle classes (to the detriment of those in the poorest and richest classes) can be detected at the asymptotic 
equilibrium. 

Dependence of the asymptotic stationary distribution on differently weighted welfare measures. The profile 
of the asymptotic stationary distribution depends on the parameter 7 characterising the welfare measures. When 7 
decreases and all other data are kept unchanged, at the asymptotic equilibrium an increase of the fraction of individuals 
belonging to the middle classes can be detected (and, correspondingly, a decrease of those in the poorest and richest 
classes). 


In this connection we evaluated for the models under consideration the Gini index G and the tax revenue TR, whose 
definition is recalled next. Let the Lorenz curve represent (on the y axis) the cumulative percentage of the total income 
of a population earned by the bottom percentage of individuals (on the x axis). The Gini index G corresponds to the 
ratio A 1 /A 2 of the area A 1 between the line of perfect equality y = x and the Lorenz curve and the total area A 2 under 
the line of perfect equality. It takes values in [0,1], 0 representing complete equality and 1 maximal inequality. The 
tax revenue is the total amount of tax collected in the unit time and redistributed as welfare provisions: 


n n n— 1 

tr = y, y, y pm s%k 

h= 1 k= 1 j= 1 


w i x i 


(E"=t w i X i 


X h x k , 


( 16 ) 


with Xj denoting the fraction of individuals in the ;-th class at equilibrium. 


4 Social mobility and its correlation with economic inequality 


This paper focuses upon an important socio-economic issue: 


social mobility]^] 


4 


The expressions social mobility and economic mobility are used here to designate the same concept. 




4.1 Introductive remarks 


Social mobility gives a measure of the chances individuals have to improve their economic status through education 
and work. It is known that a high level of inequality in a society is especially harmful when it is complemented by 
a diminution of social mobility, see e.g. in El. Social mobility has been historically an important component of the 
“American Dream” (and of its emulations), thus making the high inequality level present in the US more tolerable. 
Today, several observers point out a diminution of social mobility, as seen for instance from the growing difficulties 
that the sons and daughters of the poor and middle classes have to access the best universities and the best paid 
jobs. As mentioned in the introduction, quantitative empirical surveys show that mobility is negatively correlated with 
inequality, i.e. more unequal societies also tend to have less mobility. 

Studying the correlation between economic inequality and mobility is of particular interest because statistical eval¬ 
uations of the economic mobility of a society are more difficult than measurements of the income distribution. This 
is due to the fact that they require to follow the evolution of the individuals’ income for at least one or two genera¬ 
tions. Correspondingly, in micro-to-macro theoretical models of economic exchanges based on kinetic equations, the 
income distribution depends only on the asymptotic equilibrium solution of the equations, while mobility estimates 
also involve the detailed structure of the transition probabilities of the model. 

In our models there are no generations and it is impossible to compare the income of one individual with that of 
his or her parents. Individuals are supposed to live and work indefinitely, or at least for the time it takes for the system 
to reach its dynamical equilibrium. During this time the income distribution evolves from some initial conditions to 
a stationary state, in correspondence to which the number of individuals belonging to a certain class remains constant 
in time: at equilibrium, the total rate of individuals leaving that class is equal in absolute value to the total rate of 
individuals arriving from other classes. This means that in the equations ([7]), whose right hand side can now be written 
in a simplified form in view of ( fl5j ), one has X,- = 0. Equivalently, the equilibrium populations are solutions of the 
nonlinear equations 

I I (c hk + TL ] (X))x ll X k -X i = 0 for i £ {1, ...,n}, 

h=lk= 1 

with C hk and T' hk (X) as in |8j) and ( jToj ). 

The equilibrium populations can be obtained numerically as asymptotic solutions of the full evolution equations. 
We have evidence from several numerical solutions that the approach to equilibrium is exponential in time, in the 
sense that the “Cauchy convergence norm” converges exponentially to zero as t — >■ °° (see El). 

Below we introduce quantitative indicators suitable to define the mobility in our models. Then, we check the 
correlation of mobility with inequality. To this end, we analyze the dependence of mobility on parameters such as 
the taxation rates gap and welfare means-testing; we find interesting relations between the mobility and the inequality 
indicators. 


4.2 Mobility in our model: detailed calculation and dependence on the model parameters 

The mobility evaluated in this paper is always meant at equilibrium, when the rates of the individuals entering and 
leaving an income class are the same (but of course differ from class to class). We do not consider long-term trends, 
characterized by an increase or diminution of the total income of the society. In our models the social mobility is 
essentially a short-term mobility, which gives each individual a probability to go up or down the income ladder. A 
large value of the mobility should be regarded as positive for the economy, because it acts as a stimulus and as a reward 
for the individuals who are more capable and hard-working. While in the average the upward mobility and downward 
mobility are the same at equilibrium, we can imagine that those who go up are economically stronger than those who 
go down. This might be evidenced, in an improved version of the model, by introducing some heterogeneity factors, 
and so labelling the individuals not only with their income class i but possibly also with some “fitness” index related 
to education, number of economic links, etc. 

To start with, we define the probabilities for an individual belonging to the income class i, for i = 2, 1, to be 

promoted to the upper class following a direct interaction or an indirect interaction (i.e., through welfare provisions) 
respectively as 

S « 

Pi,exchanges(individual) ^ Pk.i (^ T) (17) 

n+i - n k=l 



and 


.wel fare (ind ivid ual) 


E 52 Ph,k^kX k X k . 


(18) 


n+i-n z? J=1 W jX jh=lk=l 

We restrict our attention to all classes but the first and the last one, in order to avoid possible boundary effects. By 
summing the probabilities (171 and (18) we obtain the total probability [^(individual) f° r a single individual advancement 

^(individual) .e.xchan gesiind ivid ual) T 7 \,welfare(individual)i 


which is the base for the calculations that follow. 


We then define “averaged” class probabilities of being promoted to the upper class (the one-step transition proba¬ 
bilities for the class i) respectively due to direct or indirect interactions as 


(19) 


and 


for i = 2 , ...,n — 1. 


' i,welfare(class) 


1 

(l-!i-!„) 


S w/Xi 

n+i - n E" = i wjXj 


n 


E 


n 

E Ph.k^kXhXk 

k= 1 


( 20 ) 


We have first introduced similar quantities in our paper D3, where we computed the ratio between the term due to 
binary exchanges (the quadratic part of the evolution equation) and the term due to tax redistribution (the cubic part). 
This ratio gives a measure of the “liberalism” of the society. We found that with equally distributed welfare the ratio 
stays, quite surprisingly, almost constant when we pass from the middle classes to the rich classes, although one might 
expect welfare provisions to be irrelevant for the class advancement of the super-rich. Let us now call 


^i(dass) ^i , exchanges ( class) T ^i.we!fare(dass) - 


the sum of the probabilities (191 and ( 201 . 


If we make an histogram of Pi( m dividual) > n a typical situation with Gini index G = 0.368, we see that it increases 
with income and attains its largest values in correspondence of the richest classes (see Fig. 0. We know, however, 
that in real world income distributions, the fifth quintile generally includes a consistent share of the classes, since the 
richest classes are scarcely populated and most of the population concentrates in the low- and middle-income classes. 
This also happens in our models provided the initial conditions are such that most of the wealth initially belongs to the 
low and middle classes. For instance, in the example with the parameters of Fig. |T] the histogram of mobility is shown 
in Fig. [2] So, while some aspects of the mobility of the rich are interesting in principle, because they are related to the 
formation of the Pareto fat tail in the distribution, the value of the mobility of the rich classes does not have a strong 
influence on the average mobility of all classes. 


Our main measure of mobility expresses the collective probability of class advancement of all classes from the 
2-th to the (n — l)-th one. Indeed, we consider 


n— 1 n 


^exchanges(collective) / 1 \ Pk.i (1 

(1 — Ai — X n ) r,+i — Ti -_2 £ =1 


( 21 ) 


and 


and take the sum 


^welfare(collective) 


S Wi Xi 


(1 X\ X„) r i+ 1 r t i=2 YJj=\ w jXj h=i k= x 


E E Ph,k^kXhX k 


XI Pexchcmges(collective) T ^*welfare(collective 

as the parameter which represents the mobility. 


( 22 ) 

(23) 


Now, we change slightly one of the parameters of the model, in order to increase equality and decrease the Gini 
index. For instance, we pass from taxation rates 30% — 45% to 30% — 60%, leaving unchanged the parameters defining 
the direct exchanges and the welfare redistribution parameter (here we have 7 = 0.5, which means uniform welfare 
distribution). This leads to a Gini coefficient G = 0.338, to a variation A Pj( c i ass ) as in Fig. |3| and to an average variation 
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Figure 1: Individual probability of single-class advancement Pni n di V idual) obtained from a numerical simulation with 15 income 
classes. The probability has been computed from the equilibrium population and from the coefficients p^, T;, according to Eq. 
|7|, in a case where T m ;„ = 30%, T max = 45% and the resulting Gini index is G = 0.368. Note that the rich classes have the largest 
individual mobility. The initial conditions for the numerical simulation are such that most of the wealth initially belongs to the low 
and middle classes. 



Figure 2: Histogram of the class mobility Pj( c i ass ) with the same parameters as in Fig. 1. This is obtained from the histogram of the 
individual mobility, weighing each column with the normalized population of the corresponding class. Since the low and middle 
classes are much more populated than the rich classes (the first 7 or 8 classes typically comprise 4/5 of the entire population), the 
weights affect strongly the relative contributions of the individual mobilities of Fig. 1. The total mobility M is determined by a 
weighted average. 


AM = 5.7 • 10 5 . It is clear from the histogram that the mobility variation for the middle classes is positive, and that it 
prevails in the average. The histogram of the individual mobility variations shows an even sharper diminution of the 
mobility of the rich classes, see Fig. [4] In fact, one could verify that also the percentage mobility variation is always 
larger for the rich classes, both for P(jndividual) and P( c iass)- A similar behavior is found when the negative variation 
of the Gini index is obtained with a change of the welfare redistribution parameter, as to make it more favorable to 
the low and middle classes. Several trials of this kind confirm the positive correlation between equality and mobility. 
In the next section we shall consider variations of the Gini index and mobility in correspondence to contemporary 
variations in tax rates and welfare distribution parameters, which again confirm this positive correlation. 


4.3 Gini index G and mobility M in dependence on two variables: level lines and correlation 

If we compute the values of G in dependence on the two variables At = Z max — T„„„ and y, we can identify on a 
cartesian plane At — y a family of “level lines of G” corresponding to constant values of G. Similarly, we can find 
in the same plane the level lines of M. In principle, one would expect the two sets of level lines to be distinct. As a 
matter of fact, it is impossible to write an explicit functional dependence of G and M on At and y, and the definitions 
of the two quantities are completely different. In particular, G depends only on the X t , while M also depends on 
the encounter frequency rates /?/,*., which are largely arbitrary. What is remarkable, however, is that the level lines 
of G and M coincide up to ca. 1 part in 10 3 , at least in the region of the plane At — y which we have explored 
through the simulations. The increments of G and M are always opposite, so the negative correlation between them 
is strongly confirmed also in simultaneous dependence on these two variables. Fig. [5] represents, for instance, the 
three coincident level lines of G and M which start from the points (AT,y) equal to (0.15,0.20), (0.15,0.15) and 
(0.20,0.15). The corresponding values of At, y, G and M are reported in Tab. EH3 























































Figure 3: Variation A//( c / ajJ ) of the class mobility when the tax rates change from 30% — 45% to 30% — 60%, and the Gini index 
correspondingly decreases from G = 0.368 to G = 0.338. 





Figure 4: Variation ^(individual) °f the individual mobility when the tax rates change from 30% — 45% to 30% — 60%, and the 
Gini index correspondingly decreases from G = 0.368 to G = 0.338. 


^min ^max 

At 

r 

G 

M 

30-45 

0.15 

0.20 

0.341 

0.002700 

25-50 

0.25 

0.28 

0.341 

0.002704 

20-55 

0.35 

0.36 

0.341 

0.002707 

Table 1: Data for level line A. 

^min ^max 

At 

r 

G 

M 

30-45 

0.15 

0.15 

0.338 

0.002712 

25-50 

0.25 

0.24 

0.338 

0.002714 

20-55 

0.35 

0.32 

0.338 

0.002717 


Table 2: Data for level line B. 


^imin ^max 

At 

r 

G 

M 

27.5 - 47.5 

0.20 

0.15 

0.335 

0.002723 

25-50 

0.25 

0.20 

0.335 

0.002723 

20-55 

0.35 

0.28 

0.335 

0.002727 


Table 3: Data for level line C. 


4.4 Mobility in the kinetic model of a physical gas 

Finally, it is interesting to compare the economic mobility defined and computed in the previous sections with the 
corresponding intuitive notion applied to a gas. We can base this comparison on the concept of temperature. Indepen¬ 
dently from its exact definition, mobility should be related to the probability that molecules pass to a state of higher 
energy after interaction with other molecules. It seems therefore reasonable to assume that mobility increases with 
absolute temperature. In order to establish a correlation with inequality, we should then determine how the “Gini index 
of a gas” varies with temperature. For the Boltzmann-Gibbs exponential distribution, the Gini index is independent 
from the temperature and equal to 0.5. A particular attention is deserved by the K'-generalizcd distribution introduced 
by Kaniadakis ED. It takes into account relativistic kinematics and is therefore suitable to describe hot gases and 
plasmas, and at the same time it provides an excellent fit for income distributions both of real world ED and of kinetic 
models m. The Kanidakis distribution has a Gini index G that increases with temperature. The analytical formula 
which gives G as a function of the parameters a and K of the distribution contains four Gamma-functions and is given 
in El. When a = 1 and K = 0, the distribution reduces to an exponential function with G = 0.5. For any other value 
of a, G is an increasing function of K (compare the examples of Fig. [6). When the Kaniadakis distribution is applied 

































to a relativistic gas, K is related to the gas parameters by 


1 . me 

= 1 + 


knT 


(24) 


where m is the atomic or molecular mass, c the speed of light, k\> the Boltzmann constant and T the absolute temper¬ 
ature EM- Eq. ( |24{ i shows that K increases when T increases. It follows that a higher temperature leads to higher G. 
This amounts, for the gas, to a negative correlation between equality and mobility, which disagrees with the results of 
empirical economic studies and of our kinetic model. 


This comparison points to one of the many qualitative differences between physical and economical interactions. 
In physical systems it happens quite often that larger gradients in density, pressure, etc. cause a quicker dynamical 
evolution, so that, in some sense, “greater inequality causes greater mobility”. In economy this “stimulus effect” has 
been theorized to work for certain growth phases, see |jT9l and references therein, but is apparently not the dominant 
factor at equilibrium. Instead, mobility is positively related to equality, at least as an average over all income classes. 
We have seen that the mobility of the low and middle classes is the key element under this respect. 



Figure 5: Level lines in the plane At — y corresponding to the following values of Gini index and mobility. A: 
G = 0.341, M = 0.00270. B: G = 0.338, M = 0.00271. C: G = 0.335, M = 0.00272. Numerical data for the three 
lines are given resp. in Table |TJ [2j [3j 

G 



Figure 6: Dependence of the Gini index of the Kaniadakis distribution on the parameters a and K, for K in the range 
[0,1] and the four values a = 1,2,3,4. 


5 Conclusions 


In this paper microscopic models of economic exchanges, taxation and redistribution, described by kinetic equations 
are considered. Some features of the emerging income distribution curves are discussed. In particular, we compare 










models with fiscal systems characterized by a different spread between the minimum and the maximum tax rate and 
models in which the welfare provision is weighted differently for different income classes. The main contribution of 
the paper is represented, in our opinion, by the introduction of a measure of social mobility. This is possible due to 
the structure of the equations, which involve several transition probabilities for the passage from an income class to 
others. Together with the Gini index, this measure serves as a tool to investigate the correlation between economic 
inequality and social mobility. The correlation turns out to be clearly negative: a lower mobility corresponds to a 
higher inequality. This is in agreement with what is observed in real world situations. While this negative correlation 
has been the object of several studies and debates in the socio-economic arena, to our knowledge it has not been treated 
within the brand of mathematical and physical approaches to complex systems. 
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